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ABSTRACT 

We study the emission from the hot interstellar medium in a sample of nearby late 
type galaxies defined in Paper I. Our sample covers a broad range of star formation 
rates, from ~ 0.1 M Q /yr to ~ 17 M /yr and stellar masses, from ~ 3 • 10 8 Mq to 
~ 6 • 10 10 Mq. We take special care of systematic effects and contamination from 
bright and faint compact sources. We find that in all galaxies at least one optically 
thin thermal emission component is present in the unresolved emission, with the av- 
erage temperature of (kT) — 0.24 keV. In about ~ 1/3 of galaxies, a second, higher 
temperature component is required, with the (kT) = 0.71 keV. Although statistically 
significant variations in temperature between galaxies are present, we did not find any 
meaningful trends with the stellar mass or star formation rate of the host galaxy. The 
apparent luminosity of the diffuse emission in the 0.5-2 keV band linearly correlates 
with the star formation rate with the scale factor of Lx/SFR s» 8.3 • 10 38 erg/s per 
M /yr, of which in average ~ 30 — 40% is likely produced by faint compact sources 
of various types. We attempt to estimate the bolometric luminosity of the gas and 
and obtained results differing by an order of magnitude, log(Lboi/SFR) ~ 39 — 40, 
depending on whether intrinsic absorption in star-forming galaxies was allowed or not. 
Our theoretically most accurate, but in practice the most model dependent result for 
the intrinsic bolometric luminosity of ISM is Lboi/SFR ~ 1.5 • 10 40 erg/s per M /yr. 
Assuming that core collapse supernovae are the main source of energy, it implies that 

esN ~ 5 • 10~ 2 (-Esat/10 51 ) of mechanical energy of supernovae is converted into 
thermal energy of ISM. 
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1 INTRODUCTION 

It is well established that the number of high-mass X- 
ray binaries and their collective luminosity scale with 
the star formation rate (SFR) o f the host galaxy 
llGrimm. Gilfanov. fc Sunvaevl |2003l ; iLehmer et all l20ld : 
iMineo. Gilfanov. fc SunvaevlbOlC Paper I hereafter). This 
fact is well understood in terms of the short evolutionary 
time scales of high- mass X-ray binaries (HMXBs), both the- 
oretically (e.g. Verbunt fc van den Heuvel 19 95) and obser- 
vationally (e.g. IShtvkovskiv fc Gilfanovl [20071 ). 

Star-forming galaxies are also known to possess sig- 
nificant amounts of hot ionized gas of ~sub-keV temper- 
atures, which is a source of copious X-ray emission (e.g. 
iGrimes etT al. 2005; Tiill mann et al.ll2006l ; lOwen fc Warwickl 
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120091 ; iLi fc Wand |2012h . The morphology of the d iffuse 
X-ray emission sugg ests (e.g. IStrickland et all l200Ch that 
the gas is in the state of outflow, driven by the collec- 
tive effect of supernovae and winds from ma ssive stars (e.g. 
IChevalier fc Cleggj|l985l h iTvler et al.l l|2004l ) demonstrated 
that the spatial distribution of soft diffuse X-ray emission 
correlates with the sites of recent star formation in spiral 
arms as traced by mid-infrared and H a emission. The to- 
tal gas luminosity was found to generally corr elate with the 
star f ormation rate of the host galaxy (e.g. | Grimes et al.l 
120051 ; I Owen fc Warwickl [20091 ; ILi fc Wangj|2012l ). Although 
the majority (if not all) of the studies agree on the presence 
of such correlation, its slope and the scale factors determined 
by different authors differ considerably. One part of this dis- 
crepancy could in principle be caused by the difference in 
the galaxy samples (e.g. ULIRGs versus more normal spiral 
galaxies) or analyzed regions (e.g. halos versus disks). But 
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Figure 1. Count rate of diffuse emission after the source removal versus the radius of the source region used for the point source removal, 
for the three spiral galaxies chosen as representative of our sample: NGC 3310, NGC 3631 and NGC 5775. The curves in each panel are 
normalized to give the same count rate at a source cell size of 2.5. The dotted line (green) indicates the count rate in 0.5-2 keV band, 
the dash-dotted line (red) represents the count rates in 2-8 keV band, the dot solid curve (blue) shows the count rate in 0.5-8 keV band. 
The source region radius is expressed in in units of -Rt)o%pSF ( see Sect. 1 3. fl for details). Note that NGC 3310 has the smallest, by a factor 
of ~ 1.6, value of the _D25 major radius. 



at least some part of the discrepancy is caused by the differ- 
ence in the data analysis methods, in particular in the way 
the contribution of compact sources was separated from the 
truly diffuse emission. 

The goal of this paper is to study the properties of the 
hot interstellar medium (ISM) in the sample of nearby star- 
forming galaxies from Paper I, with the main emphasis on 
the relation between X-ray luminosity of gas and star for- 
mation rate of the host galaxy. We aim to notably improve 
the previous studies in terms of accuracy of the treatment 
of systematic effects (such as spill-over of counts from bright 
sources) as well as in terms of numbers of galaxies for which 
uniform Lx — SFR data is available. 

The structure of the paper is as follows. In Section 2 we 
briefly describe our sample. In Section 3 we present the tech- 
niques used to isolate the hot ISM and obtain its X-ray spec- 
tra and in Section 4 we consider the remaining contribution 
of faint discrete sources to the unresolved X-ray emission. In 
Section 5 we present the results of the spectral analysis of 
diffuse emission. The X-ray luminosity of the ISM emission 
is obtained in Section 6 and its relation with the SFR of 
the host galaxy is investigated. In Section 7 we discuss our 
results. 



2 SAMPLE SELECTION AND DATA 
PREPARATION 

The goal of the present work is to study the unresolved X- 
ray emission from star-forming galaxies. To achieve this goal 
the latter must be separated from the emission of compact 
X-ray sources. We constructed a sample of galaxies based 
on the "primary sample" defined in Sect. 2.1 of Paper I. We 
selected galaxies whose full extent lies within a single chip of 
the Chandra detector. As the background correction is nor- 
mally performed on a chip-to-chip basis, the latter criterion 
allows a straightforward and accurate background level es- 



timation and subtraction, addressing both the instrumental 
and cosmic X-ray background (CXB). We selected and an- 
alyzed only single observations having exposure-time longer 
than 20ks. The thus constructed sample includes 21 galaxies 
(Table[3]) whose star formation rates span a range from ~ 0.1 
to ~ 20 Mq yr _1 , similar to that of the primary sample of 
Paper I. 

The data preparation was done following the standard 
CIACQ threads (CIAO version 3.4; CALDB version 3.4.1), 
exactly in the same way as described in Sect. 3 of Paper I. 



3 ISOLATING THE EMISSION FROM THE 
HOT ISM 

The truly diffuse emission from the hot interstellar medium 
is polluted by several contaminating factors. These factors 
are: (i) "spill-over" counts from bright resolved compact 
sources that have been removed from the image; (ii) un- 
resolved faint compact X-ray sources; (iii) instrumental and 
cosmic X-ray backgrounds. In order to estimate the X-ray 
emission from hot ISM, the contribution of these compo- 
nents needs to be computed and removed. In this and the 
next section we describe the procedures used to isolate the 
truly diffuse emission from contaminating components, to 
obtain the ISM spectrum and calculate its luminosity. 

3.1 Compact X-ray sources and their residual 
counts 

The point sources analysis procedures used in Paper I were 
tuned to produce a sample with optimal point source de- 
tection sensitivity, controlled incompleteness and accurate 
point source photometry. In contrast to that, the point 
source handling procedures used in this study are aimed 

1 http: / /cxc. harvard. edu/ciao3.4/index. html 
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to minimize the contamination from spill-over counts from 
bright sources, to minimize contribution of unresolved emis- 
sion from faint compact sources and to facilitate its accurate 
estimation for further correction. 

With this in mind, we start with the determination of 
the optimal size of the region to be used to remove the 
point sources counts from the image. To this end we se- 
lected three spiral galaxies, representative of our sample: a 
face-on galaxy with strong diffuse emission (NGC 3310), a 
face-on galaxy with weak diffuse emission (NGC 3631) and 
a edge-on galaxy (NGC 5775). To each galaxy, the following 
procedure was applied. We used the wavelet-based source 
detection algorithm wavdetect (see Sect. 3.1 of Paper I for 
details on the parameter settings) to search for point like 
sources in the soft (0.5 - 2.0 keV), hard (2.0 - 8.0 keV) and 
total (0.5 — 8.0 keV) energy bands. For each source, we used 
the information about the shape of the point spread func- 
tion (PSF) in the source position to determine the radius 
of the region containing 90% of source counts, following the 
method described in Paper I. In the following, we refer to 
this radius as i?go%psF- 

Based on the source lists obtained in each energy band, 
we created a set of source regions having radii ranging from 
0.5 ■ i?9o%psF to 2.5 ■ i?9o%psF with a step of 0.1. A cor- 
responding set of diffuse emission images was created for 
each galaxy, in each energy band, adopting the following 
method. Using the CIAO task dmfilth, we removed the 
source regions from the image and filled in the holes left 
by the source removal with pixel values interpolated from 
surrounding background regions (POISSON method). For 
the background region we used a circle with radius 3 times 
the radius of the source region. In order to avoid biases in 
the interpolation, we ensured that the chosen background 
annuli did not contain neighboring point sources. For each 
background region listed in the input file, we subtracted all 
the overlapping neighboring point source regions and merged 
them into a single source removal region. 

For each such image we estimated the count rate within 
the D25 ellipse using the CIAO task dmextract. The depen- 
dence of the count rate on the radius of the source region is 
shown in Fig.[TJ In all three galaxies, the count rate appears 
to follow the same trend - a rather sharp decrease at small 
source radii ^ Rgo% psf is followed by a much slower trend 
and the curves become nearly flat at R > 1.5-Rgo% psf- Some 
residual decrease beyond R ~ 1.5i?go% psf is observed in the 
case of NGC 3310, which is likely caused by the fact that at 
large values of R the source region size becomes comparable 
with the characteristic angular scale of the surface bright- 
ness variations. Note that this galaxy has the smallest, by 
a factor of 1.6, value of D25 major radius among the three 
galaxies selected for this analysis. In any case, the remaining 
variations of the diffuse emission count rate do not exceed 
~ 10% level even for this galaxy. 

Based on this analysis we adopted the source region 
radius of R = 1.5i?go%psF- This value is a reasonable com- 
promise, minimizing the contamination of diffuse emission 
by point source counts without compromising the statistics 
for the diffuse emission itself. 
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Figure 2. An example of the analysis performed to minimize the 
contamination of the diffuse emission by spillover counts from 
bright point sources. The figure shows the full band (0.5 — 8 keV) 
surface brightness profile centered at NGC1313 X-l (solid line, 
red) and the radial profile of the Chandra ACIS-S PSF at this 
position (dash-dotted line, blue). The vertical dotted line repre- 
sents the default source region radius of 1.5Rgo% psf> the vertical 
dashed line is the 8" radius used for bright point sources (see Sect. 
rOlfor details). 



3.2 Spillover counts from ULXs 

The procedure to optimize the source region size described 
in the previous sub-section allowed us to minimize the con- 
tamination of diffuse emission by spillover counts from the 
majority of point sources. However, the residual counts from 
extremely bright compact sources may still pollute the un- 
resolved X-ray emission. After visual inspection, in 5 galax- 
ies (NGC 1313, NGC 4214, NGC 4490, NGC 5474 and 
NGC 7793) we detected the presence of an excess of X-ray 
counts around very bright compact sources after their re- 
moval. For these sources, the counts located outside the de- 
fault source region can make up to ~ 10% of the total count 
rate of the unresolved emission. In these cases the size of the 
source region needs to be increased. 

In order to determine the required source region size, we 
extracted the Chandra PSF in the positions of bright point 
sources using the mkpsf task. The surface brightness profile 
centered on the positions of bright sources was constructed 
for both the observed counts in 0.5-8 keV band and the PSF. 
An example of such profiles for NGC1313 X-l is shown in 
Fig. [2] As it is clear from the figure, the measured profile 
follows the PSF up to at least ~ 5 — 8 arcsec, indicating 
that the surface brightness is dominated by counts for the 
bright point source and not by the diffuse emission itself. 
Therefore, excluding a circle with radius R = 1.5 • -Rgo%psF 
for these sources is insufficient. For these sources we decided 
to use the source region radii in the ~ 8 — 10 arcsec range, 
depending on the brightness of the source. 
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Figure 3. An example of exposure-corrected background-subtracted surface brightness profiles in the soft (dash-dotted line, blue) hard 
(solid line, magenta) bands for 3 representative galaxies. The statistically significant brightness fluctuations around zero level are present 
in some case. These variations are stronger in the hard band, can have both positive and negative amplitude and do not exceed 10% of 
the background level. This level is indicated by horizontal thin lines for the soft (dash-dotted horizontal line) and hard (dashed horizontal 
line) energy bands. This level characterizes the accuracy of our background subtraction procedure. The vertical dashed line indicates the 
radius corresponding to the D25 region. 



3.3 Instrumental and cosmic X-ray background 

The instrumental and cosmic X-ray background was deter- 
mined from the background regions located on the same 
chip, outside the body of the galaxy. In order to deter- 
mine the exact locations of the background regions, for each 
galaxy we constructed exposure-corrected radial profiles, 
with no background subtraction, in the soft (0.5-2 keV) and 
hard (2-8 keV) energy bands. The surface brightness profiles 
were constructed in concentrical elliptical annuli, around the 
center of the galaxy, with shape parameters - axis ratio and 
the position angle, set equal t o the corresponding param - 
eters of the optical emission (|de Vaucouleurs et alj Il99lh . 
The profiles were extracted over the whole Chandra detec- 
tor area. All profiles showed clear flattening in the outer 
parts, indicating that no significant X-ray emission from the 
target galaxy was present there. The background region for 
each galaxy was defined as the part of the Chandra CCD 
outside the D25, 1.5 ■ D25 or 2 ■ D25 region, depending on 
the shape of the emission profile of the galaxy. 

Once the background regions were defined, we con- 
structed for each galaxy the exposure-corrected background- 
subtracted surface brightness profiles, in order to check the 
accuracy of the background determination procedure. Ex- 
amples of such profiles are shown in Fig. [3] for the same 
3 galaxies as used for the source region analysis. We found 
that in some galaxies statistically significant surface bright- 
ness variations are present at large radii, where no emission 
from the galaxy itself is expected. 

These variations are typically stronger in the hard band. 
When present, they are equally strong in the positive and 
negative amplitude, suggesting that they may be resulting 
from systematic effects. However, as the zero level was de- 
termined from the data itself, their physical origin can not 
be excluded either. Their amplitude is within 10% of the 
background level and, to be conservative, we accepted this 
as a threshold in analyzing the spatial distribution of the 



diffuse emission. For all galaxies, this threshold was reached 
near the D25 boundary of the galaxy. 



4 UNRESOLVED EMISSION AND FAINT 
COMPACT SOURCES 

4.1 Spectra of unresolved emission 

In order to obtain the unresolved emission spectra with well 
controlled contribution of unresolved faint point sources we 
defined the luminosity limit L95 corresponding to a rather 
high level of the source detection completeness - 95%. Point 
sources with luminosity above this limit were excluded us- 
ing the source regions determined in Sect. 13. H and 13^2"! The 
sources below this limit, even though statistically signifi- 
cantly detected by the wavdetect task, were left on the 
image. This method allows us to accurately estimate the 
contribution of various types of faint compact sources to the 
unresolved emission. 

The spectra of the point source-free emission were ex- 
tracted using the CIAO specextract script. To be consis- 
tent with the analysis performed in Paper I, the extrac- 
tion regions were same as defined in the Section 4 of Pa- 
per I, ranging from ~ D25 for moderately sized galaxies 
to ~ 3/4 — 1/2 D25 for the largest ones, the bulges being 
excluded when possible. We grouped the spectra in order 
to have minimum 15 counts per channel to apply the \ 2 
fitting. The nominal background spectra were extracted as 
described in the Sect. 13.31 The resulting background sub- 
tracted spectra of unresolved emission are shown in Fig. [5] 

4.2 Unresolved high-mass X-ray binaries 

The contribution of unresolved HMXBs to the X-ray dif- 
fuse emission spectrum of each galaxy was estimated in the 
statistical way. 

In order to determine the spectrum of faint HMXBs 
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Table 1. Completeness limits and luminosities of unresolved 
HMXBs. 

UNRESOLVED HMXBS a 



Unresolevd HMXB spectrum 



Galaxy 




L 9 £ 


Lq 5 


-2 keV 




L 2 _ 


-8 


keV 








(ergs 


) 


(erg 


S 


1 , 




erg s 


- 1 ) 


NGC 


0278 


2 


13 


X 


10 38 


8.35 


x 


10 37 


1 


74 


X 


10 38 


NGC 


0520 


1 


96 


X 


10 38 


7.55 


x 


10 38 


1 


57 


X 


10 39 


NGC 


1313 


8 


22 


X 


10 36 


2.11 


x 


10 37 


4 


39 


X 


10 37 


NGC 


1569 


1 


06 


X 


10 36 


1.43 


x 


10 36 


2 


98 


X 


10 36 


NGC 


2139 


1 


84 


X 


10 38 


1.16 


x 


10 39 


2 


42 


X 


10 39 


NGC 


3079 


2 


26 


X 


10 38 


7.31 


x 


10 38 


1 


52 


X 


10 39 


NGC 


3310 


1 


53 


X 


10 38 


1.50 


x 


10 39 


3 


12 


X 


10 39 


NGC 


3556 


2 


97 


X 


10 37 


2.82 


X 


10 38 


5 


86 


X 


10 38 


NGC 


3631 


1 


12 


X 


10 38 


1.03 


X 


10 39 


2 


15 


X 


10 39 


NGC 


4038/39 


1 


53 


X 


10 38 


9.55 


X 


10 38 


1 


99 


X 


10 39 


NGC 


4194 


3 


79 


X 


10 39 


5.09 


X 


10 39 


1 


06 


X 


10 40 


NGC 


4214 


1 


81 


X 


10 37 


1.98 


X 


10 37 


4 


12 


X 


10 37 


NGC 


4490 


1 


11 


X 


10 38 


6.64 


X 


10 38 


1 


38 


X 


10 39 


NGC 


4625 


1 


73 


X 


10 37 


4.54 


X 


10 36 


9 


45 


X 


10 36 


NGC 


5253 


2 


22 


X 


10 37 


4.14 


X 


10 36 


8 


62 


X 


10 36 


NGC 


5474 


1 


30 


X 


10 37 


4.34 


X 


10 37 


9 


04 


X 


10 37 


NGC 


5775 


1 


11 


X 


10 38 


1.52 


X 


10 39 


3 


16 


X 


10 39 


NGC 


7090 


8 


55 


X 


10 37 


1.44 


X 


10 38 


2 


99 


X 


10 38 


NGC 


7541 


4 


22 


X 


10 38 


1.98 


X 


10 39 


4 


12 


X 


10 39 


NGC 


7793 


3 


37 


X 


10 37 


7.99 


X 


10 36 


1 


66 


X 


10 37 


UGC 


05720 


6 


25 


X 


10 38 


6.15 


X 


10 38 


1 


28 


X 


10 39 



a Soft and hard band luminosity of unresolved HMXBs com- 
puted as described in the Section 14.21 

^ The luminosity threshold used for removal of compact 
sources. It is defined as the 0.5-8 keV luminosity correspond- 
ing to 95% source detection efficiency (see sect. 14. ill . 



we combined the spectra of all resolved point sources with 
log Lx < 37.5 detected in the galaxies of our sample. The 
chosen luminosity limit allows us to include compact sources 
that are as faint as possible without compromising the statis- 
tics of the final combined spectrum. The spectrum of each 
detected point source was extracted from a circular region 
centered on the source position given by wavdetect, the ra- 
dius of the source region was fixed at -Rgo%psF- For the 
background region we used an annulus with the inner and 
outer radius equal to Rgo% psf an d 3 • Rgo% psf respectively. 
From this annulus we subtracted all the overlapping neigh- 
boring point source regions. Average spectra of faint point 
sources in each galaxy were extracted and weighted ARFs 
and RMFs were obtained, using the tasks dmextract, mkwarf 
and mkrmf . These spectra were then further combined into 
one using the combine_spectra task. The resulting spec- 
trum, grouped to have minimum 15 counts per cannel is 
shown in Fig. [4] 

This spectrum was initially modeled with an absorbed 
power law model, however, we noted the presence of excess 
counts at low energies. We modeled the soft excess with 
two thermal components with best fit temperatures of 0.13 
keV and 0.45 keV. The other parameters of the fit were: the 
power law slope F = 1.8 and the hydrogen column density 
iV H = 4.9 • 10 21 cm -2 , the reduced chi-squared xl = 1-009 
for 125 degrees of freedom. 

In order to compute the contribution of faint HMXBs, 
for each galaxy the spectrum was re-normalized in order 



1CH 



10-' 



» 10- s 





fit 1 W w+w ~ v ** 



Energy (keV) 

Figure 4. The average spectrum of unresolved high-mass X-ray 
binaries. It was obtained by combining the spectra of all the com- 
pact sources with logLx < 37.5 detected in the resolved galaxies 
from our sample (see Sect. I4.2l for details). The best-fitting model 
includes two thermal components with temperature of 0.13 keV 
and 0.45 keV respectively and a power law with slope 1.8, ab- 
sorbed by a hydrogen column density of 7Vh = 4.9 • 10 21 cm -2 . 



to give the expected 0.5 — 8 keV band luminosity of un- 
resolved HMXBs. The latter was calculated by integrating 
the average X-ray luminosity function (XLF) for HMXBs, 
eq. (18) of Paper I, between 10 34 erg/s and the complete- 
ness limit L95 (above which the detected compact sources 
were excluded). The XLF normalization was adjusted for 
each galaxy to match the number of point sources detected 
above L95 (corrected for the expected contribution of back- 
ground AGNs). The contribution of unresolved HMXBs is 
plotted for each galaxy in Fig. [5] and presented in a tabular 
form in the Table [T] along with the values of the complete- 
ness limit, L95. 



4.3 Unresolved LMXBs, CVs and ABs 

In early type galaxies, the number and luminosity of low- 
mass X-ray binarie s (LMXBs) sca les well with the stellar 
mass of the galaxy (lGilfanovll2004l ) . 

This may suggest that the contribution of unre- 
solved LMXBs can be predicted from the stellar mass 
enclosed in the analyzed regions of galaxie s. How ever, 
as it was first noticed by Bogdan fc Gilfanovl (|2010l ) and 
M mco, Gilfanov. fc Sunvaevl (|2012l ). the scale in the Lx — 
Mi, relation for LMXBs appears to depend on the age of 
the stellar population, the Lx/M* ratio being smaller in 
younger stellar environments. Therefore the average scal- 
ing relation of (Gilfanov 2004|) will overestimate the con- 
tri bution of unresolved LMXBs . This was further confirmed 
bv lZhang. Gilfanov. fc Bogdan! (|2012l ) who found strong age 
dependence of the Lx/ A/* in 20 nearby early type galax- 
ies. However, although IZhang. Gilfanov. fc Bogdanl (I2012T ) 
demonstrated the significance of the effect, they were not 
able to separate the age dependence from the dependence 
on the globular cluster content of the galaxy, and did not 
provide age-dependent scaling relations which could be used 
to predict numbers and luminosity of unresolved LMXBs. 
On the other hand, as LMXBs are the (very) minor con- 
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AM 




NGC4194 





Figure 5. The observed spectra of diffuse X-ray emission for the sample of resolved galaxies. The models show the predicted contributions 
of unresolved HMXBs (solid curve, blue), LMXBs (dashed, red), CVs and ABs (dash-dotted, green), young stellar objects (dotted, 
magenta) estimated as described in Sect. [4] 



stituent of the compact source population in the analyzed re- 
gions, their scaling relation can not be renormalized for each 
galaxy individually, using the numbers of detected compact 
sources, in the same manner as the contribution of unre- 
solved HMXBs was computed in the Section T4.2I These fac- 
tors make the accurate calculation of the contribution of un- 
resolved LMXBs to the unresolved emission in star-forming 
galaxies currently impossible. On the other hand, one of the 
selection criteria in constructing the galaxy sample was de- 
signed to control the overall contribution of LMXBs (Section 
2.1 of Paper I). The latter was further minimized by the way 



in which the analysis regions were selected (Section 4.2 of 
Paper I). For the above reasons, similar to Paper I, we chose 
not to subtract the contribution of unresolved LMXBs. 

For the same reason, although the Lx — M* scaling rela- 
tion is well established for cataclysmic varia bles (CVs) and 
activ e binaries (ABs) in Solar ne ighborhood (Sazonov et al.l 
2006) and in early type galaxies ( |Bogdan fc Gilfanov 2 01 ll ) , 
it can not be directly applied to the stellar environment in 
star-forming galaxies. Following the same line of arguments, 
we chose not to subtract their contribution. 

In order to roughly access the amplitude of the contri- 
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bution of faint sources associated with the old stellar popu- 
lations, we plot in Fig. [5] along with the measured spectra of 
unresolved emission and with the contribution of unresolved 
HMXBs, the spectra of unresolved LMXBs, CVs and ABs, 
estimated using the average scaling relations for these types 
of sources. The LMXB spectrum was computed following 
the same method as used in the Section [4.21 for calculation 
of the contribution of unresolved HMXBs. We used the av- 
erage LMXB lu minosity function and scaling relation from 
iGilfanovl ||2004 ). corresponding to L x (> 10 37 erg/s)/M Q = 
8 x 10 39 erg/s/lO n M . For the spectrum of faint LMXBs 
we assumed a power law with photon index V = 1.8 and 
JVh = 10 21 cm -2 . In computing the combined spectrum of 
CVs and ABs we used their average 2-10 keV scaling relation 
for ellipti cal galaxies L-^/Lk = 3.1 ■ 10 27 erg/s/ Lk,q deter- 
mined bv lBogdan fc Gilfanovl (|201lh . For the effective aver- 
age spectrum of CVs and ABs we assumed a flat power law 
spectrum (r = 2), in agreement with the shape of the spec- 
tru m of unresolved emission in gas-poor elliptical galaxies 
(e.g. iBogdan fc Gilfanovl l201C0. To this spectrum we applied 
same absorption as for unresolved LMXBs, Nn = 10 21 cm~ 2 . 
We emphasize that these calculations and, correspondingly, 
the curves plotted in Fig. [5] represent an upper limit on the 



contribution of faint LMXBs. The same is possibly true for 
CVs and ABs. 

4.4 Unresolved young faint objects 

Several types of intrinsically faint X-ray objects are ex- 
pected to be found in star-forming environments. These in- 
clude protostars, young stellar objects, young stars and col- 
liding wind binaries. Although they are intrinsically X-ray 
faint, their collective radiation can, in principle, contami- 
nate the unresolved emission. Unlike ellipticals, there are 
virtually no gas-free star-forming galaxies, therefore sepa- 
ration of faint objects from the truly diffuse ISM in the 
same ma nner as it was done for fa int objects in early type 
galaxies dBogdan fc Gilfanovll201lh is difficult, if not impos- 
sible. On the other hand, these objects typically have hard 
spectra, corresponding to temperatures of the order of a 
few ke V, when approximated by br e msstrahlung emissio n 
model (|Feigelson fc Montmerleill999l ; IWinston et al.ll2007l 1. 
This is significantly larger than sub-keV temperatures, char- 
acteristic fc£_the_rSM_Jri_^ar-FOTming galaxies. With this 
in mind, Bogdan & Gilfanov (201 1]) assumed that the hard 
unresolved emission of late type galaxies is due to radia- 
tion from faint compact objects and attempted to calibrate 
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their scaling relation in the 2-10 keV energy band, obtaining 
Lx/SFR ~ 1.7- 10 38 (ergs- 1 )/(A/ Q yr- 1 ), which is about 
~ 7% of the total HMXB luminosity per SFR and com- 
parable to the contribution of unresolved LMXBs in our 
sample. We use this relation to indicate the possible con- 
tribution of faint young compact sources to the unresolved 
emission in Fig[5] along with other components. In com- 
puting it we assumed that their spectra can be described 
by bremsstrahlung e mission with temperature kT ~ 3 keV 
(I Winston et al.ll2007l ). to which we applied same absorption 
as for unresolved LMXBs, CVs and ABs, A H = 10 21 cm" 2 . 

4.5 Contamination by faint compact sources — 
Conclusions and caveats 

Out of several different types of unresolved compact sources, 
the contribution of only one, faint HMXBs, can be accu- 
rately calculated and removed from the unresolved emis- 
sion. For other types of sources (LMXBs, CVs and ABs, 
young stellar objects), observations provide some useful con- 
straints which can be used to approximately characterize 
their contribution. These predictions are summarized in Ta- 
ble [2] They are however insufficient for accurate calculation 
of the luminosity and spectrum of unresolved sources and 
their removal. 

In order to proceed further we relied on the assumption 
that the combined emission of faint compact sources of all 
types has sufficiently hard spectrum, harder than the emis- 
sion of the ISM, and that the ISM makes the dominant con- 
tribution to the strong soft component observed in the ma- 
jority of our galaxies (Fig[SJ). This assumption is justified by 
observations of many individual objects and combined spec- 
tra of gas-free early type galaxies. We thus restricted our 
analysis to the soft 0.5 — 2 keV band and did not attempt to 
draw any quantitative conclusions regarding its properties 
above 2 keV. Based on the extrapolation of scaling relations 
obtained for various types of faint compact sources in the 
hard band, we conclude that in the 0.5-2 keV band, the 
background- and HMXB-subtracted unresolved emission is 
a reasonable approximation to the emission of the diffuse 
ISM. Furthermore, along with the total passband luminosi- 
ties, we will use spectral fitting results to explicitly sepa- 
rate the thermal emission from the power law component, 
presumably associated with emission of faint sources with 
harder spectra. 

As a caveat, we note that if a diffuse gaseous com- 
ponent of a >(few) keV temperature is present in star- 
burst galaxies, it would n ot be identified in the analysis of 
iBogdan fc Gilfanovl (|201lh . hence their scaling relation for 
young compact objects would overestimate their luminos- 
ity and the presence of such hot ISM component would be 
missed in our analysis as well. Similarly, although the ma- 
jority of faint compact sources do have rather hard spectra, 
there are stellar sources characterized by sub-keV emission 
temperatures. They would be attributed to the ISM com- 
ponent in our analysis. Same is true about the supernova 
rem nants - with the te mperature in the ~ 0.5 — 1 keV range 
(e.g. lLong et all l201Ch their spectra are somewhat harder 
than the soft ISM component (see below), but they will con- 
tribute to the overall unresolved emission in the 0.5-2 keV 
range. These uncertainties appear to be currently unavoid- 
able. 



5 SPECTRAL ANALYSIS 

We analyzed the background subtracted spectra in the 0.5-8 
keV band. The spectra were approximated by a model con- 
sisting of one or, if required by the data, two thermal compo- 
nents (mekal in XSPEC) and a power law. The metallicity of 
the thermal component was fixed at the solar value. To this 
model we applied double photoelectric absorption compo- 
nent (phabs*phabs), of which one was fixed at the Galactic 
value and the second was left as a free parameter in order 
to allow for any absorption intrinsic to the target galaxy. To 
this model we added a component accounting for the contri- 
bution of unresolved HMXBs, which parameters were fixed 
at the values determined as described in the Sect. [4~2l De- 
tails of the best-fitting models for our sample galaxies are 
described in Table [3] 

In all galaxies, at least one thermal component was re- 
quired, with plasma temperature in the range of ~ 0.2 — 0.3 
keV. For about ~ 1/3 of galaxies the second thermal compo- 
nent was also statistically significantly required by the data 
(F-test probability smaller than 10 -3 ). Its temperature was 
in the range of ~ 0.5 — 0.9 keV. These spectral character- 
istics are typical of the diffuse emission from normal star- 
burst galaxies and consistent with measurements published 
by other authors for several galaxies from our sample (e.g. 
IStrickland et~aHl2004lGrimes et al.ll2005l ). 

For ~half of galaxies no intrinsic absorption was re- 
quired by the data (as before, we used the F-test probability 
of 10 -3 as the threshold). In these galaxies the intrinsic A?h 
was fixed at zero value. 

The power law component was included in order to ac- 
count for the emission from unresolved compact sources, 
other than high-mass X-ray binaries, whose contribution 
could not be accurately determined from the mass and star 
formation rate of the galaxy (see Sect(4]). For the majority 
of galaxies the power law slope was poorly constrained, we 
therefore fixed the slope to the value of T — 1.7 expected 
for low-mass X-ray binaries and let only its normalization 
be a free parameter of the fit. In four galaxies (NGC 2139, 
NGC 4194, NGC 5474 and NGC 7090), no power law compo- 
nent was required by the data. In these cases it was excluded 
from the fit. 

The luminosity of the best fit power law component in 
the soft and hard band for each galaxy is listed in the first 
two columns of the Table [2] along with the predicted lumi- 
nosities of various unresolved components calculated in the 
section [4] and their sum. All luminosities are corrected for 
the Galactic and intrinsic absorption. As it is obvious from 
the Table, the sum of the predicted luminosities of LMXBs, 
CV/ABs and YSOs often exceed the observed luminosity of 
the power law component, in agreement with the consider- 
ations presented in the Section [4.51 



6 X-RAY LUMINOSITY OF THE DIFFUSE 
ISM 

We calculated luminosity of the diffuse ISM using a few dif- 
ferent approximations, and, finally, estimated its bolometric 
luminosity. 
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Figure 6. The Lx — SFR relation for the diffuse ISM emission in the 0.5-2 keV band. In the left panel, the ISM luminosity was computed 
as the total luminosity of the unresolved emission from which contribution of backgrounds and unresolved HMXBs was subtracted (see 
section [6j. The right hand panel shows the luminosity of the mekal component(s). In both panels the luminosity is corrected for the 
Galactic absorption. 

The best fit relations, described by equation ((TJ and J2}, are shown by the solid line. The errors are given at lcr confidence level. 



6.1 The 0.5-2 keV luminosity 

First, we measured the 0.5-2 keV luminosity of the unre- 
solved emission, from which contribution of backgrounds 
and unresolved HMXBs was subtracted. For each galaxy, we 
computed count-to-erg conversion for the best fit spectral 
model, which included correction for Galactic absorption. 
Using these conversion coefficients we calculated the energy 
flux and luminosity of the diffuse emission from the observed 
count rate of the unresolved emission. Before the conversion 
we subtracted from the latter the count rate of the back- 
grounds and the predicted count rate of unresolved HMXBs. 
We measured the luminosity in the 0.5-2 keV energy band 
and named it as diffuse emission luminosity LQ'f_ 2kcV . As 
discussed in the section [4] the 2-8 keV band is likely sub- 
ject to significant uncertainties due to contribution of var- 
ious types of faint unresolved sources, we therefore do not 
attempt to perform any quantitative analysis in this energy 
band. The so obtained luminosities are listed in column 5 
of Table and plotted against the SFR of the host galaxy 
in the left hand panel of Fig. [6] The existence of a strong 
correlation is obvious from the plot. We fitted the data with 
a power law model log L x lff = log A + /3 log SFR, using the 
X 2 minimization technique and obtained the best fit slope 
P = 1.06. The error on the slope computed with the stan- 
dard Ax 2 = 1 prescription is ~ 0.01, thus, from the formal 
point of view, the best fit value of ft statistically significantly 
deviates from unity. However, as the \ 2 is very large in the 
minimum, \ 2 ~ 6500 for 19 d.o.f., the formal statistical error 
is not applicable. More appropriate is the value, computed 
via least square covariance matrix, which takes into account 



the scatter of the data points. The so computed error is 
a/3 ~ 0.1. With this, more realistic error, the slope is com- 
patible with unity. We therefore fit the Lx — SFR data with 
the linear formula Lx = A x SFR using \ 2 minimization 
technique and obtained the following relation: 

Lo i f-2kcv(ergs- 1 ) w (8.3 ±0.1) x 10 38 SFR (M e yr" 1 ) (1) 

with a dispersion of a — 0.34 dex. The luminosity in the 
above relation is corrected for Galactic absorption. The er- 
ror is la for one parameter of interest and accounts for sta- 
tistical uncertainties only. The advantage of this method is 
that it has minimal dependence on the details of the spectral 
model. 

As another approximation to the ISM luminosity we 
compute the 0.5-2 keV luminosity of the mekal component, 
corrected for the Galactic absorption (sum of luminosities, 
when the second thermal component was required by the 
data). The so obtained luminosities are listed in column 6 
of Table and plotted against the SFR of the host galaxy 
in the right-hand panel of Fig. [6] The best-fit slope for the 
Lmokai — SFR relation is ft = 0.94. As the \ 2 is still rather 
large in the minimum (x 2 = 504.4 for 19 dof), we rely on 
the error computed from the covariance matrix, op ~ 0.1 
and conclude that the slope is consistent with unity. Fixing 
the slope at unity, we obtain the following linear relation: 

L^f-Wergs- 1 ) a (5.2 ± 0.2) x 10 38 SFR (M yr' 1 ) (2) 

with a dispersion of a = 0.34 dex. In computing the \ 2 we 
approximately took into account the asymmetry of the lu- 
minosity errors by using the positive or negative uncertainty 
depending on the sign of the deviation of the data from the 
model. 



© 2012 RAS, MNRAS 000, [[]-?? 



10 S. Mineo, M. Gilfanov and R. Sunyaev 



Comparing equations ([T]) and ((2J we can estimate that 
the average contribution of sources with hard spectra to the 
0.5-2 keV luminosity of the unresolved emission is ~ 37%. 
However, values for individual galaxies show considerable 
scatter with the rms ~ 22%, i.e. about the half of this value. 
A fraction of this hard emission is associated with various 
types of unresolved faint compact sources, but some fraction 
can, in principle, be due to the contribution of the hotter gas, 
unaccounted by the low temperature mekal components. 

Finally, we note that the result of the fitting of the L x — 
SFR relation depends not only on whether the least square 
of \ 2 minimization is used, but also on whether the fitting 
is performed in linear or logarithmic space. For example, 
the least squares fit to the data in the log space, logLx = 
log A + log SFR, gives close but not same value of the scale 
factor in eq. ©, 7.1 ■ 10 38 erg/s per M /yr (5.7 ■ 10 38 in eq. 
([2])). The difference is larger for the \ 2 minimization in the 
logarithmic space (16.2 • 10 38 and 11.6 • 10 38 respectively). 
This should be taken into account in comparing results of 
different analyses. 

6.2 Bolometric luminosity of gas 

In many tasks (e.g. energy balance of the ISM or calculation 
of the supernova feedback) the knowledge of the bolometric 
luminosity of ISM is critical. However its calculation is im- 
peded by a number of factors, all related to the smallness 
of the gas temperature. Below we attempt to compute the 
scaling relation for the bolometric luminosity of ISM. We do 
it in several steps, with each step producing better approxi- 
mation to intrinsic bolometric luminosity of gas, but at the 
same time obtaining progressively more model dependent 
result. 

For the XSPEC's mekal model with temperature in the 
range kT — 0.24 — 0.71, the bolometric correction for the 
0.5-2 keV luminosity ranges from m 3 to « 1.4. For the two 
temperature spectrum with the emission measures equal to 
the average emission measures from the Table [3l the bolo- 
metric correction factor is ~ 2. Applying this correction to 
the eq.© we obtain the scale factor for the apparent bolo- 
metric luminosity: Lboi/SFR ~ 10 39 erg/s per M /yr. This 
relation is corrected for the Galactic absorption but does not 
account for possible intrinsic absorption in galaxies. 

For almost a half of galaxies a model with intrinsic 
absorption gave a statistically significant improvement of 
the fit (Table ©. The obtained values of iV H are in the 
~ (1 — 4) ■ 10 21 cm -2 range, i.e. in the range of values ex- 
pected for star-forming galaxies. For the majority of remain- 
ing galaxies, the upper limits for the intrinsic absorption are 
of the same order. It is not clear, whether intrinsic absorp- 
tion is real or is an artifact of an inadequate spectral model. 
Even if it is real and measured with reasonable accuracy in 
our spectral fits, the calculation of the intrinsic luminosity 
of the gas involves large absorption correction factors, due 
to small temperature and relatively large Nh- The final re- 
sult is sensitive to the details of the spectral model at low 
energies, below the low energy threshold of Chandra detec- 
tors. It is further complicated by the degeneracy between 
the model parameters resulting in rather large uncertainties 
in the computed luminosity. 

With these caveats in mind, we attempt to derive the 
scaling relation for the intrinsic bolometric luminosity of 
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Figure 7. The L x — SFR relation for the 0.3-10 keV luminos- 
ity of the thermal component(s), corrected for both the Galactic 
and intrinsic absorption. The best fit relation, described by equa- 
tion are shown by the solid line. The errors are given at Icr 
confidence level. 



the ISM. For galaxies for which intrinsic absorption was re- 
quired by the data, we compute the 0.3-10 keV luminosity 
of thermal component(s) corrected for both Galactic and in- 
trinsic absorption. The so obtained luminosities are listed in 
column 7 of Tableland plotted in Fig. [7] The x 2 fit to the 
data with the linear relation gives: 

^atlokcvtergs- 1 ) = (7.3±1.3) x 10 39 SFR(M yr" 1 ) (3) 

As before, we took into account the asymmetry of the error 
bars for the luminosity, which in this case notably affected 
the result, due to strongly asymmetric errors for some of 
the galaxies (Fig[7]). In computing the best fit, we excluded 
NGC 1569 (the leftmost point in Fig. [7}. This galaxy has 
the largest value of L X /SFR and the smallest value of the 
negative uncertainty, significantly smaller than other galax- 
ies in our sample, and strongly affects the best fit, shifting 
it in the upward direction by a factor of 1.5. 

Applying the bolometric correction for the average spec- 
trum, Lboi/io.3-io keV ~ 2, calculated above, we estimate 
the scaling relation for the intrinsic bolometric luminosity 
of gas: Lboi/SFR ~ 1.5 • 10 40 erg/s per M /yr. Compar- 
ing this with the apparent bolometric luminosity, derived in 
the beginning of this subsection, we conclude that correction 
for intrinsic absorption introduces a factor of ~ 10. As dis- 
cussed above, due to several unavoidable uncertainties these 
numbers should be treated with certain caution. 



7 DISCUSSION 

7.1 Comparison with previous results 

Based on yiMM- Newton observations of a s ample of nearby 
late-type galaxies, I Owen Sz Warwick] (|2009T l studied the un- 
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Figure 8. Gas temperature of the thermal components for the X-ray diffuse emission spectra plotted against the SFR (left) and the 
stellar mass (right) of the galaxy. Filled circles (grey) and squares show temperatures of the cooler and hotter component respectively. 
Note that the latter was not required in all galaxies. The errors are given at la confidence level. 



resolved emission in the 0.3-1 keV band. They measured the 
soft Lx/SFR ratio in the range from 2.5 x 10 38 to ~ 10 39 
(erg s~ 1 )/(Mq yr" 1 ). In their analysis, they removed the 
contribution of the most luminous point sources only, there- 
fore the unresolved emission may include some contributions 
of X-ray binaries. Their SFR determination was based on 
FUV luminosities and corrected for dust attenuation effects. 
For our sample of resolved galaxies there is a close correspon- 
dence between the attenuation corrected UV-based SFR and 
the total SFR estimator adopted in the present work (see 
Fig. 12 of Paper I). Therefore, to compare our results with 
those of lOwen fc Warwick! (|2009f ). we only need to convert 
the X-ray luminosities from 0.3—1.0 to 0.5—2.0 keV band. To 
do this conversion, we assumed a two component spectrum 
consisting of the optically thin plasma component with tem- 
perature kT = 0.24 keV and a power law component with 
the photon index V = 2, with the absorption column density 
of ./Vh = 3 x 10 21 cm" 2 . The ratio between the normaliza- 
tions of the power law and thermal component was fixed at 
1.36 x 10" 2 , the average value determined from the spec- 
tral fits. This led to the conversion factor of 0.86. With this 
value of the conversion factor, their Lx/SFR ratios lie in 
the range from 2.9 ■ 10 38 to ~ 1.2 • 10 39 (erg s _1 )/(M© yr" 1 ), 
in reasonable agreement with the scale factors given by our 
eq s. CD and j2l) . Furthermo re, to be closer to the derivation 
of lOwen fc Warwick! (|2009T ), we did not subtract contribu- 
tion of unresolved HMXBs and obtained a scaling factor of 
7.5 • 10 38 (ergs _1 )/(M0 yr _1 ), in yet better agreement with 
their range of Lx/SFR values. 

iGrimes et all ((2001 ) computed the Lo.3-2kov/LFiR ra- 
tios for a sample of ultra-luminous infrared starburst and 
dwarf starburst galaxies and derived a mean Lq.3_2 keV /Lfir 



ratio of 10 -4 with a dispersion of ~ 0.4 dex. Their energy 
band, 0.3 — 2.0 keV, is almost identical to the 0.5 — 2.0 keV 
band used in this work. In order to compare our results with 
the derivations from these authors, we only converted the 
FIR (42.5 - 122.5 (im) luminosity to total IR (8 - 1000 fim) 
band using a conversion factor of 1.7. The latter was 
then converted into SFR using th e value of Lir /SF R = 
2.2 ■ lO 43 (ergs- 1 )/(M yr- 1 ) from iKennicuttl l|l99ct ). We 
thus found L X /SFR » 1.3 ■ 10 39 (erg s _1 )/(M Q yr _1 ). This 
value exceeds by a factor of ~ 2 our eqs CO and ([2]) and 
is also near the high e nd of the Lx/SFR determined by 
lOwen fc Warwick! {2009). On the other hand, it is ~ 5 times 
smaller than our estimate of the abso rption corrected 0.3- 
10 keV luminosity, eq.©. We note that lGrimes et alj (120051 ) 
did not do any removal of point sources, they rather relied 
on the spectral separation of the contributions of gas and 
compact sources. 



iLi fc Wand l|2012i ) recently computed the relation be- 
tween the X-ray luminosity of the galactic coronae and 
the SFR for a large sample of nearby edge-on galax- 
ies. For the absorption corrected 0.5-2 keV luminosity of 
the thermal emission, they found Lx/SFR ~ 1.3 ■ 10 39 
(ergs )/(M0yr -1 ), larger than our eq.CO) and ([2]), but 
smaller than our estimate of the 0.3-10 keV luminosity, 
eq.([3j)- In estimating the X-ray luminosity, they removed the 
detected point sources and accounted for the contribution of 
unresolved CVs and ABs in the spectral fits. 

Thus, there is a scatter of a factor of up to ~ 2 — 5 be- 
tween the scale factors obtained by different authors. This is 
not surprising, given the difference in the calculation meth- 
ods, bright compact source removal and absorption correc- 
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tion or lack of thereof, in the galactic samples and sizes of 
the analyzed regions between different authors. 



luminosity of ISM, uncorrected for possible intrinsic absorp- 
tion in star- forming galaxies, whereas the eq.Q relies on our 
estimate of the intrinsic bolometric luminosity of the gas. 



7.2 Physical parameters of the gas 



We plot in the Fig[8]the best fit gas temperatures vs SFR 
and the stellar mass of the galaxy. As it is obvious from the 
plot, there are statistically significant variations in the gas 
temperature between galaxies. These variations have rms — 
0.05 keV and rms = 0.11 keV for the low-temperature and 
high-temperature plasma respectively. However, we did not 
find any obvious trends with the SFR and stellar mass of 
the galaxy. 

Knowing the emission measure and the volume occupied 
by gas we can estimate the gas density. Assuming that gas 
occupies a disk like volume with thickness of the order of the 
exponential s cale height of the ISM determined from its X- 
ray emission (|Li fc Wandl2012T ) we find the gas densities in 
the ~ 10" 3 - 10~ 2 cm 3 range. The radiative cooling times 
of the gas are in the ~ 10 r — 10 9 years range. The high end 
of this range, ~ 10 9 yrs, found for some of the galaxies, is 
likely to be caused by too simplified assumptions about the 
shape of the volume occupied by the hot plasma. 



7.3 Efficiency of the ISM heating by supernovae 

The existence of the linear relation between the ISM lumi- 
nosity and the star formation rate allows one to estimate 
the efficiency of the ISM heating by the supernovae, assum- 
ing that they are the main source of energy powering the 
ISM e mission (e.g.lchevalier fc Clegdl 19851; Istrickland et all 
|2000| . 12004 iGrimes et al.ll2005l ; iBogdan fc Gilfanovl |200S| ). 
Indeed, for the majority of galaxies the gas cooling times 
are relatively short, suggesting that an external energy in- 
put is required in order to maintain the X-ray luminos- 
ity of the gas throughout the star formation event. Con- 
sidering only core collapse supernovae, their rate can be 
related to the SFR: ~ 1 SN per 100 yrs per 1 M©/yr 
(|Botticella et alj l2012h . Assuming further that one super- 
nova releases Esn = 10 51 erg of mechanical energy, the rate 
of the mechanical energy release due to core-collapse super- 
novae is: Esn ~ 3.2 • 10 41 x SFR erg/s. Comparing it with 
the relation for the intrinsic bolometric luminosity of the 
gas given in the end of the Section [6] we constrain the effi- 
ciency of the conversion of the mechanical energy of the core 
collapse supernovae in to thermal energy of ISM: 



-2 / Esn \ 

V 10 51 J 



(4) 



Of course all caveat regarding the accuracy of the estimation 
of the bolometric luminosity of ISM given in the Section 
6, are fully applicable here. For example, using the scaling 
relation for the apparent bolometric luminosity, Lx ~ 1 • 
10 39 x SFR, one would obtain an order of magnitude smaller 
number, csn ~ 3 ■ 10~ 3 . 

Eq. Q di f fers b y a factor of ~few from recent results 
bv lLi fc Wand" l|2012t ). who found the energy conversion ef- 
ficiency of core collapse supernovae in the range 2.0 x 10~ 3 
to 1.2 x 10~ 2 . We be lieve that the mai n source of discrep- 
ancy is the fact that iLi fc Wand (|2012h used the apparent 



7.4 Contribution of the ISM emission to the total 
X-ray luminosity due to star formation 

Comparing the relations eq. (jTJ)- (J3j) with the results of Paper 
I we conclude that above 0.5 keV, the apparent luminosity 
of hot ISM is about « 20% of the HMXB luminosity. The 
behavior of intrinsic luminosities depends on the details of 
the spectrum of the collective emission of HMXBs at low 
energies and significance of the intrinsic absorption. With 
present numbers, assuming that there is no significant soft 
component in the HMXB spectra, the intrinsic X-ray lumi- 
nosity of the ISM may be comparable or even exceed that 
of HMXBs. 



8 CONCLUSIONS 

For the sample of nearby star-forming galaxies from Paper 
I, we studied characteristics of the diffuse emission. We took 
special care of various systematic effect and we believe that 
we are in good control of the contamination of unresolved 
emission by bright compact sources of all types and by un- 
resolved faint high-mass X-ray binaries. The contribution 
of other types of faint sources (LMXBs, CVs and ABs and 
young stellar objects) can not be accurately subtracted given 
the present knowledge of their scaling relations. We argued 
that their contribution in the soft band is not dominant. 

At least one soft thermal component is statistically sig- 
nificantly required in all galaxies in our sample. Its has the 
mean temperature of {kT} ~ 0.24 keV with rms = 0.05 
keV. In about ~ 1/3 of galaxies, the second thermal compo- 
nent is required by the data, with the mean temperature of 
{kT) « 0.71 keV and rms = 0.11 keV. Although we observe 
statistically significant differences in temperature between 
galaxies, no clear trends with the stellar mass or star forma- 
tion rate could be found. 

The luminosity of diffuse emission in the 0.5 — 2 keV 
band, as well as luminosity of its thermal component(s), 
correlate well with the star formation rate. The correspond- 
ing scaling relations are described by eqs.([T])-(|3]). We at- 
tempt to estimate the bolometric luminosity of the gas 
and obtained results differing by an order of magnitude, 
Lx/SFR ~ 10 39 - 10 40 (ergs- 1 )/(M Q yr- 1 ), depending on 
whether the intrinsic absorption was allowed for. Our the- 
oretically most accurate, but in practice the most model 
dependent result for the intrinsic bolometric luminosity is 
L bo i/SFR~ 1.5-lO 4O (ergs~ 1 )/(M yr- 1 ). This number ex- 
ceeds by about an or der of magnitude the scaling relations 
obtained previously ( Grimes et alj 120051 : 1 Owen fc Warwick] 
2009; ILi fc Wandl20ia r Various caveats associated with its 
derivation are discussed in the section 16.21 Using this scal- 
ing relation, we estimate the efficiency of supernova feedback 
into thermal energy of ISM, egjv ~ 5%. 
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Table 2. The observed luminosity of the power law component of the unresolved emission and predicted contributions of various types of faint sources in soft and hard ht^ 
bands. 
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a Soft and hard band luminosities of the power law component determined from the spectral fitting of the spectrum of the unresolved emission from which backgrounds 
and contribution of unresolved HMXBs (Table [TJ was subtracted (sect.O, corrected for both Galactic and intrinsic absorptions. 
" Predicted luminosities of unresolved CVs and ABs (Sect. l4~3l . 

Predicted luminosities of unresolved LMXBs (Sect. l4~3ll . 
^ Predicted luminosities of unresolved YSOs (Sect. l4~4t . 
e Sum of the predicted luminosities of LMXBs, CV/ABs and YSOs. 
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Table 3. X-ray spectral analysis of the diffuse emission. 
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NGC 7541 


14.7 


4.9 
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NGC 7793 
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c 



a Star formation rate from Spitzer and GALEX data (see Sect. 6 Paper I). 

" Stellar mass from 2MASS da ta (see Sect. 5 Paper I). ^ 

Inclination angles taken from Tully (1988). ^ 
^ Luminosity of the diffuse emission in the 0.5-2 keV band corrected for Galactic absorption, backgrounds and contribution of unresolved HMXBs subtracted (Sect.[6]|. 

e Luminosity of the thermal component of unresolved emission in the 0.5 — 2keV band, corrected for Galactic absorption (Sect. [BJ. ~|' 

f Luminosity of the thermal component of unresolved emission in the 0.3 — lOkeV band corrected for Galactic and intrinsic absorption (Sect. [6j|. cn 

9 Intrinsic hydrogen column density, errors are quoted at 90% confidence level. S 
^ Low-temperature plasma component of the spectral model, errors are quoted at 90% confidence level. 

1 The emission measure of the first thermal component, errors are quoted at 90% confidence level. 2 

3 High- temperature plasma component of the spectral model, errors are quoted at 90% confidence level. g 

^ The emission measure of the second thermal component, errors are quoted at 90% confidence level. ^3 
' X 2 divided by the number of degrees of freedom for the best fit spectral model. 

5 1 

^ . 

I — 1 
OX 



